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We present preliminary results for the topological charge and susceptibility determined from the low-lying 
eigenmodes of the Wilson-Dirac operator. These modes have been computed on dynamical configurations with 
Nt = 2 non-perturbatively improved Wilson fermions. We compare our results with the eigenmodes of fermions 
in the quenched approximation. 
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1. INTRODUCTION 

For several reasons there is particular interest 
in the low-lying eigenmodes of the Dirac opera- 
tor in QCD. From phcnomcnological models it is 
expected that for sufficiently light quark masses 
physics, e.g. hadron correlators, is dominated by 
these modes. Calculations on the lattice enable 
us to study the relevance of low- lying eigenmodes 
for various observables directly. 

Furthermore, the eigenmodes of the Dirac oper- 
ator carry information about the topological con- 
tent of the background gauge field. The topologi- 
cal properties can also be probed using gluonic 
methods. These however typically suffer from 
large fluctuations on very short scales and there- 
fore require some sort of smoothing procedures. 
Usually the quantities of interest are affected by 
the applied filtering method. 

Results from fermionic methods on the other 
hand are expected to be sensitive to the effects 
caused by explicitly breaking chiral symmetry 
when using Wilson fermions. Although the sym- 
metry is expected to be restored in the contin- 
uum limit, these effects are potentially large on 
coarse lattices. However, these can be reduced by 
applying Symanzik's improvement program. Re- 
moving discretization effects of 0(a), where a is 
the lattice spacing, is achieved by adding a single 
counterterm to the action. The improved Wilson- 



Dirac operator then reads: 



M = 1 — K 



(1) 



where H is the hopping term of the usual Wilson 
fermion action. 

For the calculations presented here we used 
0(20-50) dynamical configurations with Nf = 2 
flavours of degenerate quarks for six different 
combinations of f3 and K sca . We used lattices of 
size 16 3 x 32 with a lattice spacing a that varies 
between 0.104(1) and 0.095(1) fin (using the force 
scale tq — 0.5 fm). For the lightest and heaviest 
sea quarks the ratio mps/ffiv was 0.60 and 0.83, 
respectively. For comparison we performed the 
same calculations using quenched configurations 
at (3 = 6.0 and a = 0.093(1) fm. 

2. COMPUTING LOW-LYING EIGEN- 
VALUES 

The Wilson-Dirac operator M is a non- 
hermitian operator. The chirality Ui of an eigen- 
mode Xi can be calculated using the correspond- 
ing right eigenvector vf. 



i(t) = ^y^v\{x,t)l 5 Vi(x,t). 



(2) 



*Talk given by D. Pleiter at the Lattice Conference 2001, 
Berlin, Germany. 



We can distinguish between two types of eigen- 
modes. Firstly, complex modes with non- 
vanishing imaginary part and u>i = and, sec- 
ondly, real modes with vanishing imaginary part 
and Ui 7^ 0. 
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Figure 1. 0(100) smallest eigenvalues of the 
masslcss Dirac operator calculated on 20 gauge 
dynamical (left) and quenched (right) configura- 
tions. 
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Figure 2. Density of the real modes as a func- 
tion of the real part. The area under the curve is 
normalized to one. The graphs are ordered from 
heaviest, i.e. infinite, (top) to smallest (bottom) 
sea quark mass. 



In principle, it is possible to find the real modes 
by determining the zero modes of the hermitian 
operator 75M. If A = 1 — no is a real mode of 
M, then 75M will have a zero mode for k — I /a. 
In practice it might however be difficult to de- 
tect several real modes which lie close to each 
other. We therefore decided to calculate the 70- 
100 smallest eigenvalues of the non-hermitian op- 
erator M using the Arnoldi algorithm (lj . Results 
for the massless Dirac operator 

1 1/1 1\ 
D = —M - 3 

2k 2 \k c kJ w 

are plotted in Fig. [l]. 

3. REAL MODES 

For dynamical fermions one would expect 
eigenmodes around zero to be suppressed. As can 
be seen from Fig. this seems indeed to be the 
case. This effect is more pronounced if one looks 
at the density of the real modes (Fig. |J). 

Real modes close to (or below) zero, which cor- 
respond to zero modes of 75 M at k < K c , are 
believed to be the origin of so-called exceptional 
configurations, which are encountered frequently 
in the quenched approximation when approaching 



the chiral regime. Although these modes are ex- 
pected to be suppressed for dynamical fermions, 
there is a non-zero probability for this kind of 
configurations to occur (see Fig. ||). 

The correlation of the chirality Ci(r) = 
(1/V) J2 X t tyuj^x + r, t) shown in Fig. || was 
found to be narrow, confirming the observations 
made in Typically the chirality of small real 
modes is concentrated in a few time slices. 

4. TOPOLOGICAL SUSCEPTIBILITY 

The topological susceptibility normalised by 
the volume is defined by 



X 



v 



(4) 



The topological charge Qtop can be obtained via 
the Atiyah-Singer index theorem, which is ap- 
proximately true on the lattice. It establishes a 
relation between the difference in the number of 
real modes with negative (n + ) and positive (n~) 
chirality and the topological charge: 



Q 



top 



(5) 



For a subset of configurations we compared the 
results for Qtop with those obtained from gluonic 




Figure 3. Chirality per timeslice u>i(t) for real 
modes (upper picture) and the normalized cor- 
relation function of the chirality (lower picture) 
Ci(r) on a configuration at /? = 5.29, K soa = 
0.f3450. The real mode with the smallest eigen- 
value is plotted with a solid line. 

operators using the Boulder smoothing technique. 
The maximum difference found was ±2. For a 
more detailed comparison see ||. 

The results for the topological susceptibility are 
plotted in Fig. ||. The discrepancy between the 
results from the fcrmionic method applied in this 
work and from the gluonic method used by Hart 
and Teper is likely to be explained by a lack of 
statistics. Another source of uncertainty arises 
from the fact that we only calculated the lowest 
70- fOO eigenvalues. Although the probability for 
real modes is expected to decrease when moving 
away from zero, it does not become zero. To check 
for this possible error, we looked at the topologi- 
cal susceptibility as a function of the largest real 
mode taken into account (see Fig. ||). For all data 
sets we found a reasonable plateau. 

5. CONCLUSIONS 

Comparing real modes of quenched and dy- 
namical configurations we found evidence for 
an expected depopulation of the region around 
zero. The results for Qtop obtained by using the 
Atiya-Singer-theorem were in reasonable agree- 
ment with numbers obtained from gluonic mcth- 



Figure 4. Topological susceptibility as a func- 
tion of the pseudoscalar meson mass. The solid 
symbols show the results from this work, while 
UKQCD results using gluonic methods j| are 
plotted with open symbols. 
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Figure 5. Topological susceptibility versus the 
cut-off of used real modes. 

ods. Our preliminary results for the topological 
susceptibility were found to be similar to other 
calculations using the same configurations but 
gluonic methods. 

ACKNOWLEDGEMENTS 

The numerical calculations were performed on 
the Hitachi SR8000 at LRZ (Munich), the Cray 
T3E at EPCC (Edinburgh) and ZIB (Berlin). We 
wish to thank all institutions for their support. 

REFERENCES 

f. R.B. Lehoucq, D.C. Sorensen, C. Yang, 
"ARPACK Users Guide", 1997; K. Masch- 
hoff, D. Sorensen, 1996. 

2. T. DeGrand, A. Hasenfratz, Phys. Rev. D 64 
(2001) 034512. 

3. A. Hart, this conference. 

4. A. Hart, M. Teper lhcp-lat/0108022[ |. 



